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Abstract 

Layzer's approximation method for investigation of two fluid interface structures associated 
with Rayleigh Taylor instability for arbitrary Atwood number is extended with the inclusion 
of second harmonic mode leaving out the zeroth harmonic one. The modification makes the 
fluid velocities vanish at infinity and leads to avoidance of the need to make the unphysical 
assumption of the existence of a time dependent source at infinity. 

The present analysis shows that for an initial interface perturbation with curvature exceeding 
l/(2\/]4), where A is the Atwood number there occurs an almost free fall of the spike with 
continuously increasing sharpening as it falls. The curvature at the tip of the spike also increases 
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with Atwood number. Certain initial condition may also result in occurrence of finite time 
singularity as found in case of conformal mapping technique used earlier. However bubble 
growth rate is not appreciably affected. 
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Hydro dynamic instabilities such as Rayleigh Taylor instability(RTI) which sets in when a lighter 
fluid supports a heavier fluid against gravity or Richtmyer Meshkov instability(RMI) which is ini- 
tiated when a shock passes an interface between two fluids with different acoustic impedances are 
of increasing importance in a wide range of physical phenomena starting from inertial conflnement 
fusion (ICF) to astrophysical ones like supernova explosions. In ICF, the capsule shell undergoes 
the RTI both in the acceleration and deceleration phases. RTI can retard the formation of the hot 
spot by the cold RTI spike of capsule shell resulting in the destruction of the ignition hot spot or 
autoignition[T] — [1]. The hydrodynamic instabilities lead to development of heavy fluid 'spikes' pen- 
etrating into the lighter fluid and 'bubbles' of lighter fluid rising through the heavier fluid. Different 
approaches have been used for the study of such problems. Among these Layzer's[5] approach ap- 
plied to single mode potential flow model[6] — [11] is a useful one giving approximate estimate of both 
Rayleigh Taylor and Richtmyer Meshkov instability evolution. The bubbles were shown by Zhang[7j 
to rise at a rate tending asymptotically to a terminally constant velocity while spikes were shown 
to descend with a constant acceleration. However, whether for bubbles or for the spikes, Zhang's 
analysis was applicable only for Atwood number A = 1, i.e, only for fluid- vacuum interface. An 
extension to arbitrary value of Atwood number A was done by Goncharov[8j. Within limitations of 
Layzer's model as pointed out by Mikaelian|r2j bubbles were shown to rise with a velocity tending to 
an asymptotic value dependent on A and having a fairly close agreement with the simulation results 
of Ramaprabhu et al. [13] . But the spikes were found to descend with a terminal constant velocity in 
contrast to a constant acceleration as obtained by Zhang [7] for A = 1. 

Asymptotic spike evolution in Rayleigh Taylor instability behaving almost as a free fall was 
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obtained by Clavin and Williams [H] and also by Duchemin et al. [12] by conformal mapping method. 
Associated with the free fall of the spike, the surface curvature of the spike was also found to increase 
with time (i.e the spike sharpens as it falls). 

In the single mode Layzer model with generalization [S] for arbitrary Atwood number the equation 
to the interface taken in X-Y plane as 

y = ri{x,t)=r]o{t) + ri2{t)x^] (1) 

with Tjoit) > 0, ri2(t) < for bubble while ?7o(^) < 0, ri2{t) > for spike. The velocity potential 
describing the motion of the heavier fluid (density ph) and the lighter fluid (density pi) are ( gravity 
g is along the negative y direction) 

(phix, y, t) = a{t) cos (A;a;)e-^(^-''"W); y>0 {heavier fluid) (2) 

(j)i{x, y, t) = bo{t)y + bi{t) cos (fcx)e^(^-''«(*^); y<0 {lighter fluid) (3) 

where k is the wave number and a{t), bQ{t), bi{t) are amplitudes. This conventional single mode Layzer 
model has the drawback that rather than conforming to the physical requirement: Viy ^ as y ^ —oo 
it necessitates the assumption of a time dependent source at y ^ — oo[IS]. To avoid this difficulty 
we modify the single mode Layzer model by replacing the zeroth mode term bo{t)y in Eq.(3) by a 
second harmonic term viz, 

M^^y^t) = &i(t)cos(A;x)e'=(^-*W) +62(t)cos(2A;x)e=^'=(^-*(*)); y<0 (4) 
Eqs.(2) and (4) give 

Vhy{y -> oo) = -^jj^^oo = (5) 



viy{y -oo) = --^]y^_oo = (6) 

The kinematic boundary conditions at the interface (1) are 
dri dri 

^+V,.-=V,y (7) 

Q^i^hx - Vi^) = Vhy - Vly (8) 

Setting the pressure boundary condition ph = Pi in Bernoulh's equation for the heavier and hghter 
fluids leads to[7]|H][n][n] - [22] 

Ph[-^ + \{y<pH? + gri] - Pil-^ + IWi)' + 9V] = hit) - Mt) (9) 

Following the usual procedure — [22] ,i-e, expanding ri{x, t) and the velocity potentials in powers of 
{kx) and equating coefficients of {kxY, (r = 0,2) we obtain from Eqs.(7)-(9) the evolution equation 
for the RT bubbles/spikes (nondimensionalized) tip elevation = krjQ, curvature .^2 = V2/k and 
velocity ^3 = k'^a/y/kg 

(10) 

f^-(3«..i)6 (U) 

dr 12(6-i)K|-ig)] * ' 

where 



(r + 4)±v/12r + 13 

2(r - 1) pi 



and 



r = t^kg (14) 

is nondimensionalized time. 

Starting from a set of initial values > 0, .^2 < and .^3 > which correspond to the description 
of temporal evolution of the tip of the bubble we arrive at the asymptotic value (r — )■ 00) 

6 ^ -\ (15) 
and 



Masymp - ^ 3^^^ ^ 3^ > [^ja™' ~\31 + A ^^^^ 

(by classical we mean the single mode Layzer approximation as used by Goncharov[8]). Two values 
coincide as A{= ^^^^) ~^ 1- The growth rate of the development of the height of the bubble tip 
is shown in Figure 1 and compared with classical value. It is seen that presence or absence of a 
source does not give rise to any qualitatively significant change in the growth rate of the bubble 
height [8] [9] [H]. 

To get spike like behavior of the perturbation of the interface we used 
< 0, 6 > and ^3 < (17) 
Corresponding to a start from an initial value 



1 r + 1 



6 



Eq.(ll) shows that ^2 increases monotonically (-^ > for all r > 0) while from Eq.(12) it follows 
that the depth of the spike tip below the surface of separation increases continuously = ^3 < 
and ^ < 0). These are shown in Figures 2(a) and 2(b) by plotting 1^2 (t) and C,3{t) as function 
of r obtained from numerical solution of Eqs.(ll) and (12) by employing fifth order Runge-Kutta- 
Fehlberg method. The initial value taken are [^2]initiai = 1-0 and [^slmiuai = —0.5 which satisfy 
condition (17) for all the following three values of r = 2{A = |), r = 5{A = |), r = 20{A = 
The value of ^2 which represents the curvature at the tip of the spike is an increasing function of 
r for every value of the Atwood number A (Figure 2(a)). Moreover for every given value of r the 
curvature ^2 increases with A. This implies that the spike continues to sharpen with time as well as 
with increasing Atwood number and is explicitly shown in Figure 3. Figure 2(b) shows that except 
very close to the starting instant the spike descends with a constant acceleration ~ —g (i.e, nearly a 
free fall). This agrees with the conclusions [7] for Atwood number A = 1. 

The time development of spiky behavior for [^2]iniUai > and [^slinuiai < is demonstrated 

in Figures4(a) and 4(b). This is shown both for increasing A with fixed r (Figure 4(a)) and with 
increasing r for given A (Figure 4(b)). But for ^ < [^2]iniUai < |\/^ with [^slinuiai < one 
encounters development of finite time singularity i.e, ^2 — ^ 00 and ^3 — t- —00 at a finite value of r. 
The possibility of the occurrence of such an eventuality at (or near) the tip of the spike is also found 
to arise when the RT instability is addressed by conformal mapping method [23] — [21] as mentioned 
by Clavin and Williams |14j. 

Finally for a trajectory starting from {^3)initiai < and < (6)mitm/ < = one 

finds that ^2 continues to increase towards ^2 = ^, i-e, the spike continues to sharpen as time 
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progresses and its speed of fall slowly decreases in magnitude. Because of the presence of singularity 
at ^2 — |(Eq.(12)), it is not possible to continue the numerical integration towards and beyond this 
point. This is shown in Figure 5 for initial values in the domain mentioned above. 
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Figure 1: (a). Variation of .^3 against A and (b). Variation of ^3 against r, where the 
black hne gives the saturation growth rate since by Eqs.(ll) and (12)and the dashed 
hne gives the classical value of growth rate of the tip of the bubble. 

11 




Figure 2: (a) Variation of ^2 against r and (b) Variation of ^3 against r with initial 
value ^1 = -0.1,^2 = 1.0,^3 = -0.5 and r =2 (dot-dashed), 5(dashed), 20(black). 
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Figure 3: Variation of ^2 against A with initial value ^1 = —0.1, ^2 = 1-0, ^3 = —0.5. 




13 





Figure 4: Shows formation of spikes for different values of (a)^ = | (dashed), 
I (dot-dashed), ^(black) for fixed r= 2.5 and (b)r = 2.1(Dashed), 2.3(dot-dashed), 
2.5(black) for fixed A = |. 
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Figure 5: (a) Variation of ^2 against r and (b) Variation of ^3 against r with initial 
value ^1 = -0.1,^2 = 0.1,^3 = -0.1 and r =2 (dot-dashed), 5(dashed), 20(black). 
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